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The Neel magnetization of 2 + 1 D antiferromagnets is composed of quark-like spin 1/2 constituents, 
the spinous, as follows from the CP^ mapping. These quark spinous are confined in both the Neel 
ordered phase and quantum paramagnetic phases. The confinement in the quantum paramagnetic 
phase is understood as arising from quantum tunneling events, instantons or hedgehog monopole 
events. In the present article, we study the approach to the quantum critical point, where the 
quantum paramagnetic phase ceases to exist. We find that irrespective of the intrinsic spin of 
the antiferromagnet, instanton events disappear at the deconfined critical point because instanton 
tunelling becomes infinitely costly and have zero probability at the quantum critical point. Berry 
phase terms relevant to the paramagnetic phase vanish at the quantum critical point, but make the 
confinement length scale diverge more strongly for half-integer spins, next strongest for odd integer 
spins, and weakest for even integer spins. There is an emergent photon at the deconfined critical 
point, but the "semimetallic" nature of critical spinous screens such photon making it irrelevant to 
long distance physics and the deconfined spinous are strictly free particles. A unique prediction of 
having critical free spinous is an anomalous exponent rj for the susceptibility exactly equal to one. 
Experimentally measurable response functions are calculated from the deconfined spinon criticality. 

PACS numbers: 75.10.-b,75.40.Cx,,75.40.Gb,75.40.-s 



I. INTRODUCTION 

Shortly after the dawning days of renormalization 
group studies^ of thermodynamic critical phenomena 
(continuous finite temperature phase transitions), this 
work was generalized to quantum critical phenomena 
(continuous zero temperature phase transitions^) in- 
duced by tuning parameters of the underlying Hamilto- 
nian rather than the temperature. Since the mid 1970s^, 
quantum phase transitions have attracted ever increasing 
theoretical and experimental activity. Quantum critical 
behavior has been obtained from quantum fluctuations of 
the order parameter— . It is then concluded that systems 
with d spatial dimensions have quantum critical points 
identical to thermal critical points m. d -\- z dimensions 
when the time direction scales as z space dimensions. In 
this traditional approach, the quantum transition is stud- 
ied via the Wilson renormalization group in which fluc- 
tuations of the order parameter are taken properly into 
account. This is the Landau-Ginzburg- Wilson (LGW) 
approach. On the other hand, some measurements can be 
interpreted as casting doubt on such a picture^. In par- 
ticular, critical exponents are coming out different than 
what is predicted. The exponents are not those of the 
classical d -\- z theory with order parameter fluctuations 
only. 

There have been recent suggestions^ that there will be 
quantum critical physics which do not follow from LGW 
order parameter fluctuations alone^^^. The new physics 
consists of the existence low energy elementary excita- 
tions intrinsic to, and existing only at the critical point, 
which will contribute and can modify the quantum crit- 
ical properties. It was postulated that these excitations 
will be fractionalized^^^. That quantum critical points 



will have unique eigenstates is generally true as long as 
the critical propagator has an anomalous exponent. Such 
excitations are expected to be fractionalized, but they 
need not be so in all cases. These critical degrees of 
freedom provide critical fluctuations beyond those of the 
order parameter fluctuations which are usually included 
in the standard Ginzburg-Landau- Wilson (LGW) phase 
transition lore. 

There are aspects of continuous phase transitions uni- 
versal to both classical thermodynamic criticality and 
quantum mechanical criticality. Both types of transi- 
tions are characterized by a diverging length scale as it 
is impossible for a macroscopic system to qualitatively 
change behavior unless there are arbitrarily large scale 
fluctuations or correlations, either thermal, quantum or 
botbii^ii. This diverging length scale makes the critical 
properties universal and independent of microscopic de- 
tails, except for the most general details like symmetry 
and dimensionality. The diverging correlations make the 
system respond to external stimuli in a scale invariant 
manner. 

The scale invariance universal to both thermal and 
quantum transitions is characterized by critical expo- 
nents. To be somewhat more explicit, we concentrate 
in relativist ic quantum critical points, but we empha- 
size that this physics can take place in other systems. 
For such a system, which we take to be an antiferromag- 
net, we are interested in the Neel magnetization Green's 
function, or staggered magnetic susceptibility. We will 
think of a transition between a Neel ordered (antiferro- 
magnetic) and disordered (paramagnetic) phase. 

In the ordered phase the transverse Green's function or 
susceptibility corresponds to spin wave propagation and 
it has a nonanalyticity in the form of a pole corresponding 
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FIG. 1: Density of states in the Neel ordered phase. 



to such propagation: 
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Here Z{uj,k) is between and 1, and the incoherent back- 
ground Gincoh vanishes at long wavelengths and small fre- 
quencies. The pole structure of the Green's function is 
clearly illustrated when one plots the imaginary part of 
the Neel magnetization propagator as shown in figure [H 
The fact that the Green's function has a pole means that 
transverse Goldstone spin waves are low energy eigen- 
states of the antiferromagnet. At criticality, the system 
has no Neel order and thus Goldstones cannot be elemen- 
tary excitations of the system. 

staggered Susceptibility in the Quantum Paramagnetic Phase at Fixed Wavevector 




FIG. 2: Density of states in the quantum paramagnetic phase. 

In the disordered phase the Green function or suscep- 
tibility corresponds to spin wave propagation with all 
three polarizations and it has a pole nonanalyticity cor- 
responding to such propagation: 

A(co k) 

{n{-uj, -k) ' n{uj, k)) = ' ^ + ^incoh(^, k) . 

(2) 

Here A{u;^ k) is between and 1, and the incoherent back- 
ground Gincoh vanishes at long wavelengths and small 



frequencies, A is the gap to excitations in the disor- 
dered phase. The pole structure of the Green's function 
is clearly illustrated when one plots the imaginary part 
of the Neel magnetization propagator as shown in fig- 
ure [21 That this Green's function has a pole means that 
triplet or triplon spin waves are low energy eigenstates 
of the disordered antiferromagnet. For 2+1 D antiferro- 
magnets, and in general for antiferromagnets below the 
upper critical dimension, the quasiparticle pole residue 
A vanishes as the system is tuned to the quantum crit- 
ical point^^^. At criticality, triplon excitations have no 
spectral weight and thus triplons cannot be elementary 
excitations of the system. 



staggered Susceptibility at the Quantum Critical Point at Fixed Wavevector 




FIG. 3: Density of states at the quantum critical point. 

On the other hand right at criticality the response 
function below the upper critical dimension (below which 
?7 7^ 0, while above 77 = 0) has nonanalyticities that are 
worse than poles 
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as obtained from the renormalization group studies of the 
nonlinear sigma model^^^^^^^^^^. Below the upper critical 
dimension 7^ is a nonintegral universal number for each 
dimensionality. This critical susceptibility has no pole 
structure, but has a branch cut. It sharply diverges at 
CO = ck and has an imaginary part for cj > ck. The 
nonanalytic structure of the Green's function is clearly 
illustrated when one plots the imaginary part of the Neel 
magnetization propagator as shown in figure [3l Branch 
cuts in quantum many-body or field theory represent im- 
mediate decay of the quantity whose Green function is 
being evaluated. Hence the elementary excitations or 
eigenstates of the noncritical quantum mechanical phases 
break up as soon as they are produced when the system 
is tuned to criticality: they do not have integrity. The 
complete lack of pole structure and the branch cut singu- 
larity below the upper critical dimension mean that the 
elementary excitations of the quantum mechanical phases 
away from criticality, the spin waves, cannot even be ap- 
proximate eigenstates at criticality as they are absolutely 
unstable. 
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The quantum critical point is a unique quantum me- 
chanical phase of matter, which under any small per- 
turbation becomes one of the phases it separates. It is a 
repulsive fixed point of the renormalization group. As far 
as the transition from one quantum mechanical phase to 
the other is continuous, and both phases have different 
physical properties, the critical point will have its unique 
physical properties different from the phases it separates. 
The properties of the critical point follow from the critical 
Hamiltonian H{gc) {Qc is the critical coupling constant), 
which will have a unique ground state and a collection 
of low energy eigenstates which are its elementary exci- 
tations. These low energy eigenstates are different from 
those of each of the phases as long as we are below the 
upper critical dimension. As a matter of principle, all 
quantum critical points below the upper critical dimen- 
sion will have their intrinsic elementary excitations. 

We have seen that below the upper critical dimension, 
the excitations of the stable quantum phases of the sys- 
tem become absolutely unstable and decay when the sys- 
tem is tuned to criticality. The question comes to mind 
immediately: what could they be decaying into? When 
one tries to create an elementary excitation of one of the 
phases, it will decay immediately into the elementary ex- 
citations of the critical point. The critical excitations will 
be bound states of the excitations of the stable phases 
the critical point separates. These bound states could be 
fractionalized as conjectured by Laughlin^ and Senthil, 
et. al.^ ^ but they need not be in all cases. These criti- 
cal degrees of freedom are responsible for corrections to 
the LGW phase transition canon^. The intrinsic quan- 
tum critical excitations contribute to the thermodynam- 
ical and/or physical properties of the quantum critical 
system. 

In the present work, we develop these ideas to find 
if they occur in antiferromagnets. We will consider the 
disordered paramagnetic phase of 2 + 1 D antiferromag- 
nets and the approach to the quantum critical point from 
such a phase. These are described by the nonlinear sigma 
model augmented by Berry phase terms as originally dis- 
covered by Haldanei-^^^^ and developed by other si^^^^^^*": 
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g = 2^/2/ S is the dimensionless coupling constant, a is 
the lattice constant, S is the microscopic spin with h in- 
cluded (not to be confused with the Euclidean action), 
and Sb is the Berry phase term. Their effect is nonzero 
only in the disordered phase of the Heisenberg antiferro- 
magnet as first suggested by Haldane^^ and worked out 
by Read and Sachdev^^. In particular, when the micro- 
scopic spins in the lattice are half-odd integers, the para- 
magnetic ground state has a spin Peierls bond order that 



breaks the lattice symmetry and is four-fold degenerate. 
For odd integer spins the paramagnetic ground state has 
a spin Peierls bond order that breaks the lattice symme- 
try and is two- fold degenerate. For even integer spins we 
have a valence bond solid that does not break the lattice 
symmetry. All of these quantum paramagnetic phases 
have a spin triplet gapped excitation, the triplon, and a 
spin zero gapped collective mode. 

Motivated by apparent anomalies and interesting ef- 
fects in the physics of cuprate superconductors, whose 
parent state is a Mott insulating 2 + 1 D antiferromagnet, 
Laughlin suggested^ that the Neel field would break up 
into constituents "quark" spinous at the quantum critical 
point between a Neel ordered and a quantum paramag- 
netic phase. A couple of years ago, Fisher, Sachdev and 
collaborators^ suggested that such a physics indeed oc- 
curs, but only in spin 1/2, 2 + 1 D antiferromagnets. 

The quantum paramagnetic phase of antiferromagnets 
is equivalent to a "charged" CP^ spinon field whose fic- 
titious U(l) charge couples to an emergent U{1) gauge 
field or photon generated by the fiuctuations of the CP^ 
field^^^^o. The CP^ mapping (n = z^az, with a the 
vector of sigma matrices) is obtained from the nonlinear 
sigma model description of antiferromagnets in terms of 
the Neel field n^^ . Haldane discovered that in 1 + 1 and 
2+1 D, the nonlinear sigma model needs to be augmented 
by Berry phase terms in the paramagnetic or spin disor- 
dered phasei^ii^ii^, which in 2 + 1 D leads to breaking 
of lattice symmetries for odd integers and half odd inte- 
ger spinsi^ii^. If one doubts the nonlinear sigma model 
mapping of antiferromagnets in the disordered phase. 
Read and Sachdev^^, starting from the Heisenberg model, 
showed that the disordered phase is indeed described by 
"charged" CP^ "quark" spinous z coupled to an emer- 
gent photon A^. 

In 2 + 1 D, the CP^ model has important tunnel- 
ing events2i which correspond to instanton hedgehog 
events that effectively make the U(l) gauge field com- 
pact. Moreover, if we start from the appropriate Heisen- 
berg lattice description, the gauge group is necessarily 
compact. Polyakov showed^^^^^ that compact QED con- 
fines as the Wilson loop^^ obeys an area law when in- 
stanton or monopole events are included. Polyakov's 
proof did not include matter, but it is believed^^, and 
we show below, that the presence of charged matter does 
not eliminate the tunneling instanton events as long as 
the charged matter is massive. Therefore, spinous are 
confined in the paramagnetic phase and are thus closely 
analogous to the quarks of Quantum Chromodynamics. 

Fisher, Sachdev and collaborators suggested that the 
Berry phase-induced quadrupling of instanton events in- 
trinsic to spin 1/2 antiferromagnets makes such monopole 
events irrelevant at the critical point between the param- 
agnetic and Neel ordered phases^. Therefore, spinous are 
deconfined at such a critical point for spin 1/2 antiferro- 
magnets. They further suggested that at the deconfined 
quantum critical point the critical exponent rj of the Neel 
field correlator is due to the decay of the Neel field into 
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the deconfined quark spinons. These spinons wih be the 
intrinsic excitations of the quantum critical point. 

We study here the question of deconfinement, its con- 
sequences and how it occurs. We find that when decon- 
finement of the Neel field into two spin 1/2 quarks occurs, 
the critical exponent rj is exactly equal to 1 regardless of 
the origin of deconfinement. Deconfinement will occur 
whenever instanton events vanish. We find that for all 
values of the microscopic spin, integers and half-odd in- 
tegers, instanton events vanish at the quantum critical 
point where the paramagnetic phase ceases to exist. The 
vanishing of instantons happens for two reasons. The 
masslessness of the spinons at criticality screens the in- 
stanton fields, making them irrelevant at long distances. 
Furthermore, for all values of the microscopic spins^ the 
Euclidean action of instanton events becomes infinite at 
the quantum critical point due to the masslessness of the 
spinons^-^, and thus the probability of instanton events at 
criticality is zero. Hence deconfinement occurs indepen- 
dent of the value of the microscopic spin. 

On the other hand, there is some dependence on the 
microscopic spin on the quantum paramagnetic phase as 
the lattice symmetry breaking depends on the spin value 
through Berry phase terms^^. We also found further de- 
pendence on the microscopic spin as the system is tuned 
to the quantum critical point since our analysis also has 
the consequence that the confinement length will diverge 
faster upon approach to criticality for half-odd integer 
spins, next fastest for odd integer spins, and slowest for 
even integer spins as a consequence of the Berry phase 
terms relevant to the quantum paramagnetic phase. 

For the first time, we write down the effective critical 
theory and from it, calculate experimental consequences. 
At criticality, there is still a U (l)-mediated gauge interac- 
tion between massless spinons, so in principle they might 
not be free. We find that at criticality, spinons are very 
mobile because of their masslessness. Hence they screen 
very effectively their gauge interaction so they are strictly 
free at long distances. The decay of the Neel field into 
two free particles leads to a critical exponent of 77 = 1 
exactly. This result can be used as a diagnostic of de- 
confined criticality, that is, for a free deconfined spinon 
critical point we predict a critical exponent r] exactly equal 
to one. 

The deconfined critical points studied and elucidated 
here seem to be different than the 2 + 1 D Heisenberg 
critical points. It has been suggested before that these 
two different types of critical points might occur in 2 + 1 
D^. One particular suggestion is that interactions irrele- 
vant to the Neel and quantum paramagnetic phases turn 
the Heisenberg critical point into a deconfined critical 
point and there seems to be indirect numerical evidence 
for such physics^^. Before making strong conclusions one 
must wait for experimental evidence and/or further and 
more explicit numerical evidence. 



II. CP^ MAPPING OF THE 0(3) NONLINEAR 
SIGMA MODEL 



We consider the CP^ mapping of S0{3) vectors n^^^ 
n = z^az (6) 

where 



with the restriction 
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inherited from = 1. The z's are bosonic. Notice that 
the CP^ map has 4 variables. The restriction |zp = 1 
eliminates one, leaving 3 independent variables. This 
seems to be one too much as the 0(3) nonlinear sigma 
model has only 2 independent variables since one is elimi- 
nated by the nonlinear condition. This is not so as one of 
the variables is redundant because of the gauge symme- 
try z ze^^ . The z's are spinors, i.e. spin 1/2 objects. 
The z's are the quark-like spinon constituents of the Neel 
field n. 

With this mapping the 0(3) nonlinear sigma model 
partition function becomes 



Z = J VzVz^S{\z\^ - l)e 



with Euclidean action 



S=^J d^r{\d^zf-\z^d^zf)+SB. 



and Berry phase 
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The lattice regularization is necessary to define the Berry 
phase, as the Berry phase has microscopic lattice sensi- 
tivity. We can decouple the quartic term via a Hubbard- 
Stratonovich transformation leading toMi^ 



Z = j VzVz^VA^5{\z\^ - l)e- 
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Now the gauge invariance is explicit as the kinetic term 
for the CP^ fields is built up of covariant derivatives. 
Since the action is quadratic in A^, the saddle point eval- 
uation about the minimum with respect to is exact 
and can be used to go back to the original action by 
substituting the solution to the equation of motion 
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In 2 + 1 dimensions, which is our main concern here, 
the Berry phase terms can be written in terms of the 
gauge fields as^^ 



Sb = ^iTrSCsQs 



(15) 



where S is the microscopic spin (not to be confused 
with the Euchdean action), s are points in the dual lat- 
tice around which a hedgehog or magnetic monopole of 
strength Qs is centered and ("s is 0, 1, 2 and 3 depend- 
ing on which dual lattice, W, X, F, Z, the monopole is 
centered. 
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Qs is given by 
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where the normal to the surface S is closed and F^j^ is the 
Maxwell tensor for A^. Such spinon induced gauge field 
tunneling event cannot happen unless the gauge field is 
compact. Some readers might object that the nonlinear 
sigma model from which we started should not apply 
in the paramagnetic phase. For those readers we point 
out that one can start from the Heisenberg model and 
use a Schwinger boson representation on the lattice and 
one obtains a similar action, and leads to the same CP^ 
model with the advantage that in such a derivation the 
emergent gauge field is necessarily compact. In fact, such 
a derivation was performed by Read and Sachdev^^. 

We can write the delta function in the partition func- 
tion in its integral form, so that 
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where now 



S = ^J d^r{\{d,- iAM^ + iA^(^t^ - 1)} 
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Integrating the z's we obtain an effective action for the 
yl's and the A's. This effective action has a minimum at 



A = and A = constant. Varying the effective action 
with A = and A = constant with respect to A gives a 
self consistency equation for A 



47r^ m / 

= — arctan — - 
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where A = 1/a is the inverse lattice constant and = 
2m? /{go). Since A is a mass term for the z's, the non- 
linear condition \z\^ = 1 generates a mass m for the 
spinous z fields for g > gc = 47r^. gc corresponds to 
the quantum critical point where the quantum paramag- 
netic phase dies and in all probability gives rise to the 
Neel ordered phase. 

The seems to play a passive role, as it has no dy- 
namics of its own in the bare action. This is not so 
once fluctuations and dynamical corrections are calcu- 
lated. The gauge field acquires dynamics through the 
z field fiuctuations leading to an effective action with 
a Maxwell action for the gauge fields. In fact the long 
wavelength gauge field propagator gives^^ 



(^/i( ka)Ajy[ka))i — ioop 
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(KK-S^^k^). (20) 



This term is the same as would be obtained from a term 
in the action of the form 



1 



37rm 



d-'rF, 



(21) 



Hence 3m7r/4 in two dimensions. Higher order cor- 
rections will renormalize such that = mh{g, m/A), 
where h is dimensionless. Higher order corrections also 
generate terms that vanish at long wavelengths and are 
thus dropped. Therefore, the effective action is one for 
massive spinon fields z coupled to a U{1) gauge field 



S=^Jd^r{\{d^-iA^)z 



(22) 



We then have a theory of "charged" spinous coupled to 
an emergent compact U(X) gauge field or photoni^ii^i^S. 



III. EFFECTIVE ACTION FOR THE 
QUANTUM PARAMAGNETIC PHASE OF 2 + 1 D 
ANTIFERROMAGNETS 

For the paramagnetic or disordered phase of 2 + 1 D 
antiferromagnets, g > gc-, we have seen that the spinon 
fields z acquire a gap. The quantum critical point corre- 
sponds to the spinous becoming massless at ^ = ^c- We 
also saw that the adiabatic fiuctuations of the spinous 
generate dynamics for the gauge fields. We have an 
emergent compact U{1) gauge field coupled to a com- 
plex spinor representing the spinon fields. In our case 
the gauge field is compact as follows from the lattice 
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Schwinger boson derivation of the effective CP^ theory 
for the Heisenberg model. Therefore, the quantum para- 
magnetic phase and quantum critical point of 2 + 1 D 
antiferromagnets is mapped to compact 2 + 1 D QED 
with two complex matter fields with an SU{2) internal 
symmetry. The internal symmetry is effectively the orig- 
inal invariance of the nonlinear sigma model and plays 
a passive role for the most interesting properties of the 
paramagnetic phase to be developed in this section. 

If the emergent photon did not correspond to a com- 
pact U{1) group, the paramagnetic phase would have 
massive spinon excitations and massless photons. On 
the other hand, on careful thought, this seems completely 
wrong as the paramagnetic phase is fully gapped and the 
emergent photon would be gapless. It turns out that com- 
pact QED in 2+1 D, as originally shown by Polyako v^^i^^ , 
has monopole tunneling events, usually called instan- 
tons or hedgehogs, such that the ground state of the 
theory is a monopole condensate. Such a ground state 
gives the photon a mass without breaking the gauge sym- 
metry and makes the Wilson loop obey an area rather 
than a perimeter law^^, rendering the theory confining. 
Polyakov's proof did not include matter. We show below 
that including the charged matter spinon fields, the the- 
ory still confines and the spectrum is fully massive and 
singlet with respect to the U{1) gauge field. 

Before moving to study the effects of instantons, we 
remind the reader that we do not have pure QED, but 
we have the gauge fields coupled to charge spinon fields 
which satisfy the nonlinear condition \z\^ = 1 as enforced 
by the Lagrange multiplier ^A. In order to find the effect 
of the spinous on the theory, they can be integrated out 
of the partition function explicitly to yield 



Z = J V{SX)VA^ 
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with 
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Using detM = exp [tr In M] , the partition function can 
be rewritten as 



Z = J V{SX)VA^e 
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over ^A to obtain 



detM j V{SX) exp [-5XMSX] = detM^/^ (27) 

where the operator M is the inverse of the square of the 
operator —2{d^ — iA^Y/{ga)-\-2m'^/ (ga). Hence the total 
partition function after taking care of the constraint, i.e. 
after integrating over (5 A, is 
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This is equivalent to 
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or disentangling the determinant in terms of a complex 
scalar field ^ 
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Hence antiferromagnets in 2 + 1 D, which were shown to 
be equivalent to compact QED coupled to spinor spinon 
fields z with the constraint \z\^ = 1 become equivalent to 
2 + 1 D compact QED coupled to charged complex scalar 
fields once the constraint is enforced. 

Now the complex scalar fields can be integrated out to 
yield their effect on the theory. When they are integrated 
out their main effect is to screen electromagnetic fields 
by a factor which at long distances takes the form 



1 



(31) 



where / > 0. The largest contributions to the screening 
effects of the <l>, and hence the spinous, at long distances 
can be easily calculated by computing and summing po- 
larization diagrams. To one loop order, in the Landau 
gauge, we obtain for the long wavelength photon propa- 
gator 



{A^{-k)A,{k)) ^ ^ ( 1 



If. eM 



37rm / 



kfjky 



(32) 



For the moment we only consider the 5X dependent 
part of the action and its path integration. We can ex- 
pand about the saddle point ^A = and path integrate 



Hence we see that the spinous decrease the effective 
magnetic charges and thus have a diamagnetic screen- 
ing effect which we parametrize by a screening constant 
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Cm- The most important diagrams contributing to this 
screening are the bubble diagrams, which we sum to ob- 
tain 
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Therefore in the random phase approximation (RPA), 
em = 1 + e^/(37rm). We finahy obtain that the effective 
action including the effects of spinous is given by screen- 
ing electromagnetic fields by the appropriate factor Cm 
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IV. CONFINEMENT OF SPINONS IN THE 
QUANTUM PARAMAGNETIC PHASE OF 2 + 1 D 
ANTIFERROMAGNETS 

Now that we have seen that once the effects of the 
spinous is included, the effective action for the antifer- 
romagnet is equivalent to compact QED with an appro- 
priate screening factor induced by the spinous, we need 
to take into account the nontrivial tunneling effects aris- 
ing from the compactness of the gauge group^^^^^^^^. The 
monopole tunneling events considered by Polyakov create 
electromagnetic fields according to 



Vxi7 = 0, V - H = 27rqS{r-r) (35) 



where 



(36) 



and q is necessarily an integer as follows from the com- 
pactness of the gauge field. These tunneling events cor- 
respond to certain tunneling events for the spinous via 
the mapping 



-iz^djj^z 



(37) 



In fact, such spinon tunneling events have been calculated 
by Murthy and Sachdev— and found to have Euclidean 
action given by 



Sci = 2p{q) In (A/m) 
leading to the instanton fugacity 

C e-^^' exp [-2p{q) In (A/m)] 



2p(g) 



(38) 



(39) 



This instanton action depends on the charge q of the 
tunneling event. For single instanton event p = p(l) 
0.06229609. 

From the form above we see that the instanton events 
would have a Coulomb law field 



(la [r 



Hi, / 



13 



(40) 



where are the charges of the tunneling events centered 
at in 2 + 1 D space-time. If A^^ are the classical gauge 
fields corresponding to such instanton configurations, a 
semiclassical expansion around such configurations with 
gauge fields A^^ + A^ yields the effective Maxwell action. 
Berry phases and instanton Coulomb gas action 



s 

QsQt 



2e^e 



(41) 



1/2 



where F^y is the Maxwell field of A^. In order to get 
the partition function of the antiferromagnet, we have to 
integrate over the Euclidean times of the instanton events 
and sum over the instanton events centered in the dual 
lattice. The partition function for our system is given by 



Z 



VzVz'^VA^e'^'^ Zi, 



where the Euclidean non-instanton action is 

1 



d^vF^ 



(42) 



(43) 



and 



E 



dr. 



(44) 



with Euclidean instanton action 



2^ 



1/2 



(45) 



Note that this is a Coulomb gas action augmented by 
Berry phases. 



A. Mapping of the Instanton Coulomb Gas to 
Sine-Gordon Theory Including the Effects of Berry 
Phases 



As we have just seen, the instanton contributions will 
provide constant average field fluctuations in the ground 
state, and this is described by a Coulomb gas partition 
function. This Coulomb gas can be disentangled by us- 
ing the fact that the Laplacian operator is the inverse of 
the Coulomb potential. Following closely the methods 
of Polyakov, and Read and Sachdev, which included the 
effects of Berry phases for the first timei^^— , we obtain 
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^inst ~ 



exp 



Stt 



exp ' 



j Vx exp 
j Vx exp 



f 

Jo 



dr 



'dxs 


2 


1 









''"^iTrSCs^Xsjqs 



^0 



{s,t) s 



dXs 



dr 



87r2 



/ 

Jo 



dr 



^ iXs - Xtf + ^a^ 

{s,t) s 



dr 



87r2 



Xs - Xt 



E 



dXs 
dr 



2C / ^^cos[nSC^Xs] 



(46) 



where p = 0.06229609 and we sum over instantons with 
charges q = 0, 1, —1 since higher monopole charges are 
strongly suppressed as they have very smah fugacity or 
probabihty^^^^^'^^. Hence, as found by Polyakov, we see 
that the Coulomb gas is equivalent to a Sine- Gordon the- 
ory. The only difference is that the Sine- Gordon theory 
we obtain, which is centered in the dual lattice, is frus- 
trated by the Berry phase terms since the argument of the 
cosine has the phase shift ttS^s with the values 0, ttS^ 27tS 
and 37r5' depending on whether the instanton is in the 
or Z sublattice^^^^i^ 
We first study the case of even integer spin. In such a 
case the phase shift 7rS(s in the cosine of the Sine-Gordon 



theory obtained from instanton events is a multiple of 27r 
and thus equivalent to zero. 

We now consider the case of odd integer S. In that case 
the phase shift of the cosine of the Sine-Gordon theory 
is for Xs in the W dual sublattice, and an even integer 
for Xs in the Y dual sublattice. Therefore, the W and 
Y lattice are equivalent and we shall call it Y. Similarly, 
for Xs in the X or Z sublattice the cosine in the Sine- 
Gordon have phase shifts of odd multiples of tt and thus 
are equivalent. We shall call them X. The Sine- Gordon 
action then becomes 



p2 ^2 rf^ 
^ - 9— - 



Xf - Xt 



E 



dx 



X 



dr 



where = [l67r^(ma)^^] / [e^e^a^] Defining 



E 



dr 



cos xj 



Y 



X = Xi + X2 , X = Xi - X2 
and going to the continuum limit, the Sine- Gordon action can be written as 



• (47) 



(48) 



a 



(Vxi)' + 4A\2 _ 2V2AX1 (Vx2)' + ( + (^)' - M2sinxisinx2 



(49) 



Since X2 has a mass of the order of the cutoff A = 1/a, 
all gradients and time derivatives of X2 are suppressed to 
zero. We must still minimize the action with respect to 
X2- The minimization gives 



M2 



X2 



— — sinxicosx2 ^ — ^sinxi + 



8A2 



8A2 



^\ (50) 



which, when substituted into the Sine-Gordon action 
gives the effective low energy action 



y 6 CyYI 



(^mXi) 



32A2 



cos (2x1 ) 



(51) 
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after the phase shift xi ^ Xi + '7r/2. This is also a 
Sine- Gordon model and leads to the same long distance 
instanton physics and confinement as found for even inte- 
gers. On the other hand, there is a microscopic difference 
as this model with odd spins breaks the lattice symme- 
try, leading to a two-fold degenerate ground state. Both 
models, besides the usual triplon excitations once con- 
finement is introduced, will also have a massive spin zero 
mode, X ^^e first model and Xi in the latter. 

We now move to the case of odd half integral spin. In 
this case the lattice symmetry is broken and the ground 



state is four-fold degenerate. The cosine term in the orig- 
inal Sine-Gordon theory has a phase shift of zero for 
and Xs in the dual W sublattice, a phase shift of 7r/2 
for the fields in the dual X sublattice, tt for the fields 
in the dual Y sublattice, and Sir/ 2 for the fields in the 
dual Z sublattice. Actually these phase shifts are for a 
spin 1/2. For spin 3/2 the phase shifts are interchanged 
among the four sublattices, but lead to similar physics 
by relabeling of the sublattices (basically interchanging 
X and Z) . Exactly similar four- fold degeneracies happen 
for all half-odd integer spins. The action in this case is 



r 5 5: {(xf - xrf + (xf - xD' + (xr - xff + (xr - xff} 



W 



+ 



X 



dxi 



X 



(52) 



Defining 



Xy = Xi+X2- X3 



XX = Xi - X2 + X4 
Xz = Xi - X2 - X4 



(53) 



we obtain 



S - 4^-^ 



J d?r {(Vxi)' + SA^xi + SA^x^ + 2K\l - 2Axi (Vx4) 



dxi 
dr 



9X2 
dr 



1 f dx3 



2\ dr 



1 f 9x4 



2\ dr 



[sin (xi + X2) sinxs + cos (xi X2) sinx4] 



(54) 



Since X2 , Xs ^nd X4 have masses of the order of the cutoff 
A, all gradients and time derivatives of X2, Xs ^nd X4 are 
suppressed to zero. We now minimize the action with 
respect to xs and X4 to get 



Xs = sin (xi + X2) cosxs 

M2 



X4 



-^sin(xi + X2)+0^^ 
"8A 

M2 



2 COS (Xl -X2)C0SX4 



--cos(xi-X2)+0 — 



8A2 



A6 



(55) 



Substitution of these expressions in the action yields 



32A4 



sin (2xi)sin (2x2) 



We now minimize with respect to X2 to obtain 



X2 



256A6 
256A6 



sin (2xi)cos (2x2) 



sin (2xi)+0 



A12 



Substituting into the action we finally get 



(56) 



(57) 



S - 1 



8192A6 



cos (4x1 )^ (58) 
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To summarize, we see that in general instanton effects 
lead to a Sine-Gordon action of the form 



s, 



87r2 



(59) 



where 



ns) 



for even integer S 
for odd integer S 
for half odd integer 5* 



(60) 



Finally, if we make the shift x ~^ X\/47r^/ [e'^CmfiS)]^ we 
obtain 



xf -M\S) cos ixHS))} (61) 



with 



h{S) = 



M\S) 



(62) 



fiS)-l 



We see that the effects of instantons are described by 
a Sine-Gordon theory. The conclusion that instanton ef- 
fects can be described by a Sine-Gordon theory is true for 
all microscopic spins. The only difference is in the mass 
of the Sine-Gordon theory we just wrote and in the factor 
inside the cosine of the Sine- Gordon theory. Hence the 
long distance confinement consequences of tunneling ef- 
fects are equivalent independent of the microscopic spins 
with just irrelevant numerical differences. Even though 
we did not discuss it in detail here, the Berry phases 
do make a difference for the paramagnetic ground state, 
which leads to breaking of lattice symmetries'^' They 
make the ground state quadruply degenerate for half-odd 
integer spins, doubly degenerate for odd integer spins, 
and a non-degenerate valence bond solid for even inte- 
ger spins. But as far as the instantons the long distance 
confinement physics is the same regardless of microscopic 
spins. 



B. Confinement of z Spinons in the Quantum 
Paramagnetic Phase 

In order to see that the effects of instanton fluctua- 
tions lead to conflnement, we now closely follow Polyakov 
and calculate the correlation or Green's function between 



electromagnetic flelds. The electromagnetic fleld is de- 
flned as before 



(63) 



The instantons, or the instanton charge density, acts as 
a source for these electromagnetic flelds via 



H,{x) = \j 



\x - y\- 



or in momentum space 



H^{k) = 27ri^p{k). 



(64) 



(65) 



It proves very convenient to introduce sources for our 
instantons. In particular we calculate 



S p(x)r}(x)d^x\ _ 



Z[ri{x)] 
Z[0] 



(66) 



where 



Zinst Vl] 



^qa5{x - Xa) 

a 

j Vx^^p J d^r 
- M^S) cos ixHS))]} . 



2 i^i^ (x - v)Y 



(67) 



Correlation functions of p can be derived from the par- 
tition function by taking derivatives with respect to r]: 
d^'^^ Zinst/ dri^'^\ In particular, we obtain 



{p{k)p{-k)) = e- k^{x{k)x{-k)) 



(68) 



The correlation function for the electromagnetic fields 
is given by 



{H^{k)H^{-k)) = {H^{k)H^{-k)) 



(0) 



'^{p{k)p{-k)) 



(69) 



where {H^{k)Hj,{—k))^^^ is the Green's function without 
instantons. With the appropriate screening, this Green's 
function is 



1 



1 

err 



kjj^ki/ 



A;2 



(70) 



The k = pole reflects the masslessness of the photon. 
Now we calculate the instanton-instanton density corre- 
lation function. Before doing so we note that in the dilute 
gas approximation the x coupling constant is extremely 
small {ma <C 1) and thus provides only a small renor- 
malization without changing the qualitative behavior of 
the theory. Hence the Sine-Gordon cosine potential may 



11 



be approximated to second order in x to high accuracy 
and 



{X{k)x{-k)) 



fc2 + {h{s)M{s)y 



(71) 



We find that the instanton charge density correlator is 
given by 



{p{k)p{-k)) = 



{h{S)M{S)y + fc2 

2 U2 



_ {h{s)M{s)yk 

~ + {h{S)M{S)f 
The electromagnetic propagator becomes 
{H^{k)H,{-k)) = 



1 



fc^fc. , k^k, {h{S)M{S)y 

fc2 



fc2 k-^ + {h{S)M{S)f 



(73) 



+ {h{S)M{S)y 



We see from the pole in the last term that the photon 
acquired a mass Mp = h{S)M{S) without spontaneous 
breaking of U{1) symmetry. This is a sign of confine- 
ment and that the ground state is a gauge singlet. If the 
reader is unsatisfied with this correct, but indirect con- 
clusion, we can check that the theory indeed confines by 
calculating that the Wilson loop^^^^^ 



(74) 



gives an area law^^^^^. Explicitly, using Stokes' theorem 
we have 



F[C] = (e*^^'^'^^'^) = (e^-^A^M^'^M) 
which, using equation can be written as 
F[C] = (^e^ i ^^^^P^^^^^^) 

with 

1 



r]{x) = 



\x-y[ 



(75) 



(76) 



(77) 



Using and dSTj) we find 
1 



F[C] = 



Pxexp 



m 

M\S) cos ixHS))]} . 



[d,{x-v)r (78) 



In the saddle point approximation we obtain 



F[C] - exp 



{-/ 



(79) 



M\S) cos {xciKS))]) . 



The classical field is obtained by solving the equation 

d\xci -V) = M\S)h{S) sin ixcihiS)) (80) 
which takes the form 

d\ei = 27r5Hz) eA{xy) + M\S)h{S) sin {xcih{S)) 
OAixy) = il ^^^^ . (81) 



(72) 

This has solution 



_ f4arctan(e-^(^)''(^)^) 
|_4arctan(e^(^)'^(^)^; 

It now immediately follows that 



z>0 
z <0 



■ (82) 



F[C] e 



-l(S)A 



7(5) =2] J' - - (83) 



+ / dzM'{S) cos {xci{z)h{S)) 

J —00 

where A is the area of the x?/-plane. 

We see that the Wilson loop satisfies the area law. The 
interaction energy between charges is given by 



E{R) = -f{S)R 



(84) 



and the theory confines for all values of the microscopic 
spins. The force between spinous is constant and equal to 
the "string tension" ^{S). The scale below which the the- 
ory confines 1/^(6') = M{S)h{S) is spin dependent but 
nonzero. Basically spinous are confined at length scales 
larger than ^(5) and are not relevant to the physics of 
the system. The low energy excitations are spin 1 triplons 
{z^az) and a spin collective mode corresponding to the 
Sine-Gordon field x- 



QUANTUM CRITICAL DECONFINEMENT 
OF SPINONS 



We saw in the previous section that spinous are con- 
fined and that the confinement arose from the topological 
tunneling events originating from the compactness of the 
theory. On the other hand, it has recently been suggested 
that there is a new class of quantum critical points whose 
properties are controlled by the deconfinement of spinous 
rather than fiuctuations of the order parameter^^^. In 
particular, it was proposed that this new kind of quan- 
tum critical points occur in 2 + 1 D antiferromagnets. 
Finally, it was further proposed that, because of the de- 
pendence on microscopic spins of the Berry phases, the 
deconfinement occurs only for spin 1/2 systems^. In the 
present section, we will study how deconfinement occurs 
and under which conditions. 
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The quantum critical point occurs when the spinous 
become massless, m = 0. As we saw above, the only 
dependence on the microscopic spins, as far as long dis- 
tance confinement properties, appears in the parameter 
of the Sine- Gordon theory that describes the nontrivial 



tunneling effects in 2 + 1 D antiferromagnets. As shown 
in the previous section, when m > antiferromagnets 
confine spinous for all microscopic spins. The presence 
of confinement is characterized by the nonzero value of 
the photon mass or confinement scale 



Mp{S) =M{S)h{S) = 



\ 



(85) 



such that confinement occurs for all energies smaller than 
Mp{S). It now follows that the confinement energy scale 
vanishes at the quantum critical point for all values of 
the microscopic spin. Therefore, spinous are deconfined 
at the quantum critical point independent of the micro- 
scopic spin. On the other hand, there is a dependence of 
microscopic spins on the approach to the critical point. 
For half-odd integer spin systems the confinement length 
scale i{S) = l/Mp{S) diverges faster than for odd inte- 
ger spin systems. The confinement length scale diverges 
faster for odd integer spin systems than for even inte- 
ger spin systems. Therefore the only dependence on the 
microscopic spin is how fast one reaches the deconfined 
quantum critical point, but deconfinement occurs for all 
microscopic spins. 

The conclusion that deconfinement is independent of 
microscopic spins is new and unexpected. Therefore we 
will check that this is indeed so with more care. At face 
value the only conclusion that can be reached is that at 
m = 0, M{S) = 0. Let's calculate the Wilson loop for 
the case of M{S) and see if the theory confines or 
not. As we saw in the previous section, to calculate the 
Wilson loop we need to evaluate 



F[C] -expj-y^d^ 

- M\S) cos ixcMS))]} . 

where 

d\xci -V) = M\S)h{S) sin ixcihiS)) . 
When m = 0, M{S) = and we have 

F[C]~exp|-i J <fr [3^ (x,; - 7?)f | 

with 

d\Xci - r/) = . 



(86) 



(87) 



(88) 



(89) 



Now the Wilson loop is calculated very straightforwardly 
by integrating by parts to obtain 



0/ 



F[C] exp - / {xci - V) iXci - r/) = 1 . (90) 



Since F[C] = e~^^, the long distance force between 
spinous is zero and they are deconfined. This result is 
independent of microscopic spins as M{S) for all 
microscopic spins. 

There is one last but equivalent way to see that de- 
confinement occurs independent of microscopic spin, i.e. 
that the instantons or compactness effects disappear at 
the quantum critical point. If we go back to the partition 
function from which all the confinement physics followed 



Z 

SzA 



4e2e^ J 



d^rFL 



-i^TrSCsQs 



1 Li?, 



-Sin 



(91) 



■ RtV + (t« - ftY^ 



1/2 



c 



m 
AJ 



2p(g) 



( is the instanton fugacity as calculated before^^ and 
p{q) > except when we have no instantons, where by 
definition p(0) = 0. ( = 1 for this last case. It now 
follows immediately that at the critical point C = ex- 
cept in the absence of instantons, K = O^q = 0. There- 
fore, Zinst = 1, there are no Berry phase terms and no 
monopole events and the theory is deconfined indepen- 
dent of microscopic spins. 



VI. SPINON DECONFINED CRITICAL 
THEORY 

As we obtained in the previous section, at the quantum 
critical point, when the spinon mass m vanishes, the in- 
stanton and Berry phase terms disappear from the theory 
and the quantum critical point corresponds to massless 
spinous coupled to an emergent U{1) gauge field with no 
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effects of topology. The emergent gauge group is noncom- 
pact at criticality. We had calculated before the effects of 
the spinous on the emergent photon and found that when 
the spinous were integrated out, we obtain "dielectric 
screening" of the emergent electromagnetic fields. This 
is easily understood as the spinous have a gap and hence 
the theory is a "semiconductor" as far as the emergent 
electromagnetic properties. The spinous are "semimetal- 
lic" at criticality and we need to calculate the screening 
properties in this case anew. 

The critical theory is described by the action 



2 

ga 



(92) 



As before, 6\ is a Lagrange multiplier to enforce the 
constraint \z\^ = 1. Just as the gauge fields acquire 
dynamics, spinon fiuctuations give dynamics to the La- 
grange multiplier 6\. The dynamics for 5\ follows from 
the same one loop spinon fiuctuations that gave dynamics 
to the gauge fields. To one loop order 



{5\(k)8\{-k)) 



7r2 V2Ay k 



(93) 



This term is the same as would be obtained from expand- 
ing the exponential of the action 



7r2 V2A/ 



(Fk 1 



j5\{k) 5\{-k) 



(94) 



(27r)3 k 

The 1/k means that 5\ fluctuations will be suppressed at 
long wavelengths and hence irrelevant to the low energy 
physics. For completeness, we keep these terms in the 
effective action, which leads to a 5\ propagator 



(95) 



Now we must calculate the screening effect of the mass- 
less spinous. This is easily done via the renormalization 
group by integrating high energy spinon degrees of free- 
dom and computing their effects on the renormalization 
of the electric charge. In this way we obtain the beta 
function 



d 2 



(96) 



This function is not dimensionless as usual because in 
2 + 1 dimensions the electric charge has dimension of 
square root of mass or momentum. This equation can be 
integrated easily and yields the renormalized charge at 
scale /i <C A 



(97) 



Therefore, from = e^/e^, we obtain the screening 
factor at scale /i = A: in momentum space 



(98) 



Spinous interact through photon exchange. The effec- 
tive interaction between spinous is given by the photon- 
photon Green's function, which we calculate in the 
Lorentz gauge 



V = {A^{k)A,{-k)) = ^ 



knki/ 



which in real space is formally written as 



7T 



e^{x — x'Y 



(99) 



(100) 



Therefore, the effective theory at criticality is described 
by the action 



/ 



2A 



- \ J d\d\'J,{r)V,^{r-r')Ur') 

+ j d^rr{r)z\r)(j''z{r) 
where 

J^{r)=i{z^d^z-zd^z^) . (102) 

We have added an external source J^(r) that couples to 
the Neel field = z^a^z because we want to study the 
{n^{r)n^ (r')) correlator. 

In order to analyze the critical action, we introduce 
a Hubbard- St r at onovich real vector field to write an 
action linear rather than quadratic in J^. The action for 
the nonlinear sigma model at criticality is then 



9A 

— \d^z\^^i5\{\z\^-l) 



+ \j dh<fr'B^{r)V-J{r - r')B,{r') (103) 
+ ^y" <frB^{r)J,{r) + j d^rJ%r)z^{r)a''z{r) . 



with 



kjj^kiy 



)fc2 



(104) 



Of course, in the partition function we must integrate 
over z, z^ ,5X and 5^. The Bfj^^s are gauge fields with the 
appropriate gauge invariance as can be seen from their 
action. We can now integrate the spinon fields z to obtain 
the partition function 



jvB, 



V5X e"^^" 



(105) 



where S'efr is, in momentum space for convenience 

5eff = exp [-J ^ yB^{-k)V-J{k)B^{k) - i5\ 
- TrlnM) 



(106) 
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The operator M is 



M{qi,q2) = [ifj^5\ 



9? 



) 5{q\ - q2) 



+ B^i-qi+q2) 
+ B^{qi - q2) 



qi + g2 
2 

q2 



(107) 



We can now move to calculate the Neel field correlator 
{n''{k)n\-k))5{k + q) = {n" {k)n'' {q)) 



1 



Z dJ''{k)dJ\q) 



1 



Tr((j° 



VB^V5\ 



d^qicfp 



Z" " ' J -"'^--■y (27r)6 
X M-\qi,p)M-\p-k,q + qi) 



-Seii 



(108) 



2 6" 



VB^V5\ 



cfiqi(fp 
Z J J (27r)6 

X M-\qi,p)M-\p-k,q + qi) 



j=o 



j=o 



To lowest order we find 



(n"(fc)n''(-fc)) = 25' 



cfp 



ab f _£l 

J (27r 



)3p2(p_fc)2 



— (- + —\^ —0.25 



k V8 



(109) 



We point out that this is what is obtained from non- 
interacting, free spinons. That is, the Neel field decays 
into two free spinons, which are then reabsorbed to re- 
constitute the Neel field. We see then that the Neel field 
Greenes function has anomalous exponent r] = 1. This 
is the unique consequence of decay of the Neel magne- 
tization into free spinons at criticality. In our model, 
the spinons are expected to be free as they can inter- 
act through emergent photons, but at the critical point 
the electromagnetic potentials acting between the charge 
currents generated by the spinons are screened from the 
Coulomb form to l//c, leading to a retarded inter- 
action ~ + '^^j^ ^ Coulomb law ~ 
Hence the spinons are expected to be free. Let's see if 
this is indeed so. The one loop correction coming from 
the effective spinon interaction is given by 



{n''{k)n\-k)) 



1 loop 

d^q2d^p 



p- (k^p) 



—0.104329. 

k 



(27r)6 q2{p-q2yp^{p^k)^{k^p-q2)^ 

(110) 



Therefore, including the interactions to first order we find 
that the spinons still behave as free as the anomalous 
exponent rj = 1. The only effect of the interactions was to 



renormalize the non-universal constant of proportionality 
multiplying 1/k. 

Our last result is a one photon interaction result, but a 
simple exercise in Feynman diagrams using the Feynman 
rules that follow from the effective action (|103p shows 
that higher order corrections lead to corrections propor- 
tional also to 1/k. Inclusion of Lagrange multiplier 
lines or photon lines internal to the spinon bubble 



gives contributions proportional to 1/khy direct compu- 
tation or by simple power counting. In these diagrams, 
the double solid lines correspond to Neel field n lines, 
the single solid lines corresponds to spinon z or z^ lines, 
the wiggly line corresponds to a photon propagator and 
the lightning bolt corresponds to SX lines or propagators. 
Similarly, inclusion of further internal lines to the bubble 
gives contributions proportional to 1/k. There are also 
corrections given by the diagrams 



which give contributions proportional to 1/k. Inclusion 
of higher order bubbles or corrections internal to the 
spinon bubble also yield contributions proportional to 
1/k. One nontrivial point which could invalidate the con- 
clusion that ?7 = 1 is if internal self energy corrections to 
the spinon propagator include In k terms. Explicit calcu- 
lation of the spinon self energy corrections shows that the 
log terms that are generated are multiplied by factors of 
k'^ which make them vanish as /c ^ and are thus irrel- 
evant to the universal low energy critical physics. Hence 
we see that at criticality, spinons are not only deconfined 
but behave as noninteracting at long distances and lead 
to critical exponent r] = 1 equivalent to that of decay of 
the Neel field into critical free spinons. 

The deconfined critical points studied and elucidated 
here seem to be different from the 2 + 1 D Heisenberg 
critical points. It has been suggested before that these 
two different types of critical points might occur in 2 + 1 
D^. One particular suggestion is that interactions irrele- 
vant to the Neel and quantum paramagnetic phases turn 
the Heisenberg critical point into a deconfined critical 
point and there seems to be indirect numerical evidence 
for such physics^^. In references'^, evidence for a contin- 
uous transition between a valence bond ordered param- 
agnet and its corresponding Neel ordered phase was pre- 
sented in both Heisenberg and XY systems. Evidence 
for a relatively large value of r] was also presented. It has 
been suggested that these transitions correspond to the 
deconfined critical points^ studied in the present work. 

We have shown that if a deconfined critical point exists 
in 2 + 1 D antiferromagnets, it will occur irrespective of 
the microscopic spin of the system, and that the confine- 
ment length will diverge faster upon approach to critical- 
ity for half-odd integer spins, next fastest for odd integer 
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spins, and slowest for even integer spins as a consequence 
of the Berry phase terms relevant to the quantum param- 
agnetic phase. Deconfinement occurs because instanton 
or monopole events vanish at criticality^^ and together 
with them, the Berry phase terms vanish toc^^. We also 
find that the emergent photon at criticality is screened 
strongly at long distances making the spinous behave as 
if they were strictly free at long wavelengths. Finally, 
the Neel critical correlator behavior follows from decay 
into free spinous, which universally leads to the critical 
exponent r] = 1. This is a diagnostic of deconfined criti- 
cality, that is, for a free deconfined spinon critical point 
we predict a critical exponent 77 exactly equal to one. 



VII. SOME EXPERIMENTAL CONSEQUENCES 
OF CRITICAL SPINON DECONFINEMENT 

Now that we have found the effective critical theory 
of deconfined critical points, we will briefly obtain some 
experimental consequences of the free spinon critical the- 
ory. As we have seen, the Neel critical propagator or 
susceptibility is given by 



{uaiuj, k)nb{-u;, -k)) = 
C 



SabG{u;,k) = Sab 



1 (111) 



where C is a dimensionless constant. Experimental mea- 
surable quantities are usually proportional to the density 
of states, which at finite temperatures is given by 



D{uo,k,T) 



— ImG(a;, k) — 

TT 1 



1 



g-a;/T • 

The imaginary part is easily calculated to be 



Im k) = - 

to give a density of states 
D{u,,k,T) = ^e{\co\-\k\) 



fc2 



(112) 



(113) 



o/T\ 



The dynamic structure factor is directly proportional 
to the density of states, S{uj, k) oc -D(w, fc, T). Hence the 
inelastic neutron scattering intensity is given by 



I{u), k) oc 



-e(\uj\ 

TT V 



\k\ 



p (1 



o/T\ 



For neutron energy losses which are a lot larger than the 
temperatures of the system, the dynamic structure factor 
is independent of T, with small corrections which are 
suppressed by exponentials e~^/^. For neutron energy 
losses which are small compared to the temperatures, the 
dynamic structure factor is 



x{u,k)o^^e{\u\-\k\^ 



1 



(114) 




Dynamic Structure Factor at Fixed Temperature (T = 1) and Wavevector = 5) 



Dynamic Structure Factor at Fixed Te 



rature (T = 10) and Wavevector (l< = 5) 



Dynamic Structure Factor at Fixed Temperature (T = 1 00) and Wavevector (k = 5) 





FIG. 4: Inelastic neutron scattering intensity as a function of 
frequency for fixed momentum and temperatures 1/100, 1/10, 
1, 10 and 100 respectively. 
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Product of Temperaiture and Dynamic Structure Factor 



T I<a)/T, k/T) 




Product of Temperature and Dynamic Structure Factor 



T lico/T, k/T) 




FIG. 5: Scaling of the product of the inelastic neutron scat- 
tering intensity and temperature with the ratios of frequency 
and momentum to temperature {oo/T and k/T). 



We plot the neutron intensities as a function of co at 
fixed T and \k\ in figure [H The plots can be summa- 
rized because the product of temperature and inelastic 
neutron scattering intensity, or structure factor, is a uni- 
versal scaling function of frequency over temperature and 
wavevector over temperature (i.e., it satisfies a law of cor- 
responding states), which we plot in figure [51 

Another important measurable quantity is the result 
of static neutron scattering experiments. The neutron 
intensity in this case will be proportional to the static 
structure factor. That is, it is proportional to the inelas- 
tic neutron response integrated over all frequencies 



I{k) (X- [ dcoOf 

J -A ^ 

- 1 

^ J\k\ 

_l 

^ J\k\ 



1 



1 



dcj _ = 

\k\ \Juj^-P 

duj = 

\k\ \Juj^-P 



1 1 



(1 



olT\ 



(1 



olT\ 
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i 1.6 




FIG. 6: Equal time neutron scattering intensity as . 
of the ratio of cutoff to momentum. 



function 



The frequency integrals have been cutoff to regularize 
ultraviolet divergences. The static neutron intensity is 
temperature independent and given by 



X(^) (X — In 

TT 




. (116) 



We plot the static neutron intensity as a function of the 
ratio of cutoff to wavevector transfer (A/|^|) in figure [H 

In some of these antiferromagnetic systems, some of 
the nuclei making up the material have nonzero nuclear 
spin and nonzero hyperfine coupling to the electronic 
spins that make up the antiferromagnet. In these sys- 
tems one can perform nuclear magnetic resonance (NMR) 
experiments. The NMR relaxation rate 1/Ti will be pro- 
portional to the local density of states. Therefore we 
have 

1 



1 - e-"/^ " 



(117) 



At frequencies small compared with the temperature, the 
NMR relaxation rate becomes independent of uj and lin- 
ear in temperature. In fact for |cc;| <C T 



;i (x2Tsgn(c^) . 



(118) 



For frequencies a lot larger than the temperature, the 
NMR relaxation rate is proportional to the magnitude of 
the frequency, and independent of temperature. In fact, 
for IcjI >T 



1 



(X 2\uj\ . 



(119) 



We point out that the NMR relaxation rate, or local sus- 
ceptibility, divided by the temperature is a universal scal- 
ing function oiuj/T. We plot the NMR relaxation rate as 
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Relaxation Rale at Fixed Temperature (T = 1/100) 





ition Rate at Fixed Tennperature (T = 10) 




FIG. 7: Nuclear Magnetic Resonance (NMR) Relaxation Rate 
as a function of frequency for temperatures 1/100, 1/2, 1, and 
10 respectively. 




FIG. 8: Scaling of ratio of nuclear magnetic resonance (NMR) 
relaxation rate to temperature with a;/T. 



the Neel ordered phase, it appears that using our meth- 
ods we cannot say much about the critical exponents that 
follow as the Neel phase perishes. If the critical point 
where Neel order disappears is indeed a deconfined criti- 
cal point as the ones studied in the present work, we can 
predict relations between the magnetization exponent f3 
and the correlation length exponent u. 

The Josephson correlation length in the Neel or- 
dered phase satisfies the renormalization group (KG) 



equatioi> 



with solution 



11,12,26,27,28 



dg_ 



^(m,^) -exp 



C(m,<?) =0 



dg' 



Pig') 



where 



Pig) M 



d/j, 



(120) 



(121) 



(122) 



is the usual KG beta function and gc is the critical value 
of the coupling constant that separates the Neel ordered 
and paramagnetic phases: (3{gc) = and gc > 0. Near 
the critical point but on the Neel ordered phase, the cor- 
relation length scales as 



e(M,5)-(5c-5)^/'^'(«°) 



1 



{9c - g) 



(123) 



a function of frequency for fixed temperature in figure [3 
We summarize the relaxation rate plots by graphing the 
scaling plot of the ratio of the relaxation rates to tem- 
perature, which is a universal function of a;/T, in figure 

IHl 

We have just obtained some experimentally mea- 
surable response functions that follow from having an 
anomalous dimension = 1 for the Neel magnetization 
propagator. Since the CP^ mapping of the Neel magne- 
tization into spinous proves quite cumbersome to study 



with /3'{gc) = d/3/dg\g=g^. Therefore the correlation 
length exponent is 

(124) 



/3'i9c) 



Similarly, the Neel magnetization a satisfies the KG 



equatior> 



11,12,28 



a{g) = (125) 
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FIG. 9: Feynman diagrams representing processes contribut- 
ing to the noninteracting partition function. The diagram on 
the top- left represents the emergent gauge field fluctuations. 
Lagrange multiplier fluctuations are represented by the 
diagram on the top-right, while spinon fluctuations are repre- 
sented by the diagram on the bottom. Wiggly lines represent 
gauge fields, dashed lines represent SX fields, and straight lines 
represent spinous. 



with solution 

a{g) = Mexp 



1 r (i+7(go) 

Pig') 



dg' 



(126) 



M is an arbitrary nonuniversal constant and j{g) is the 
anomalous dimension obtained from the magnetization 
renormalization factor Z via 



dlnZ 



(127) 



The critical anomalous dimension is then given by 77 = 
l^Qc)-, which is 1 for deconfined spinous. Near the critical 
point but on the Neel ordered phase the magnetization 
scales as 



<t(m,5)~(5c-5)-(^+^(^=»/'''^'(^=)1 
Therefore the magnetization exponent is 

W'{9c) ' 

which leads to the exponent relation 

(1+^ 



(128) 



(129) 



(130) 



For deconfined quantum critical points the correlation 
length exponent v and the magnetization exponent (3 sat- 
isfy the relation 




(3=v 



(131) 



FIG. 10: Feynman diagrams representing processes contribut- 
ing to the interacting partition function. These are interac- 
tions among spinous, Lagrange multiplier fields and emer- 
gent gauge fields. Straight lines represent spinous, wiggly lines 
represent emergent gauge fields, and dashed lines represent 6\ 
fields. 



which is a unique prediction for deconfined critical 
spinous. 

We can also approximate and predict the behavior of 
the specific heat at the quantum critical point. Since 
this is a finite temperature property, it is usually calcu- 
lated in Euclidean time, where the imaginary time direc- 
tion is finite and goes from to the inverse temperature 
P = 1/T. The effective action ()lQ3p is separated into in- 
teracting and noninteracting parts. The noninteracting 
part includes all quadratic terms. The rest of the terms 
are grouped in the interacting part. 

Before analysing the action, as with all gauge field the- 
ories, we must fix the gauge for in order to remove 
the gauge redundancy. This is done by the Faddeev- 
Poppov procedure, but in this case in which we will 
choose the Lorentz gauge, the Faddeev-Poppov deter- 
minant only provides an irrelevant renormalization con- 
stant. In this case, choosing the Lorentz gauge has the 
only consequence that the terms k/j^B^ in the action are 
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zero. Hence 



Physical Response to External Magnetic Field 



(^) ~^ 2 ^ 



(132) 



In order to compute the specific heat, we must com- 
pute the free energy which being the generator of con- 
nected Green's functions, it is equal to — TlnZ, where 
InZ is the log of the free field partition function plus 
the sum of connected vaccuum processes. The partition 
function receives contributions from the interacting and 
noninteracting parts of the effective action. This divides 
it in a product of a noninteracting Zq and interacting 
Zj partition function, such that Z = ZqZj. The rele- 
vant Feynman diagrams contributing to the noninteract- 
ing partition function Zq are shown in Figure [9l These 
contributions yield a term in the specific heat which is a 
positive constant times as expected for noninteracting 
relativistic particles. 

The leading contributions to the interacting partition 
function Zj are shown in Figure [TOl They give renormal- 
izations of the T^, noninteracting specific heat. These 
contributions also give a low temperature correction to 
the specific heat coming from interactions. Such correc- 
tion is proportional to a positive constant times T^lnT. 
The specific heat divided by is plotted in Figure [TTl 




|k|/T 



co/T 



FIG. 12: Physical response of the system to an external mag- 
netic field as a function of the ratios of frequency and mo- 
mentum to temperature. 



The magnetization is given by 



Specific Heat / Temperature 




d\nZ 1 dZ 



dBi 



ZdB^ 
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(134) 



OSe 



B=0 



FIG. 11: Ratio of the specific heat to 



and the susceptibility is 



Another interesting experimental signature of decon- 
fined, essentially free spinous at a quantum critical point 
is the behavior of the system in the presence of an exter- 
nal magnetic field. The magnetic field couples to the Neel 
magnetization vector and rotates it such that the sigma 
model partition function in the presence of the external 
magnetic field is 



J VnS{n^ - 1) e-^^ 

Sb^^ dr j d?x (133) 
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B=0 



The expectation value of the magnetization in the ground 
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state is zero as it consist of 
d 



/ - A i9Bdr'n{x ,t') • [B{x ,t') x n(x ,r )) 
oBi[x^T) L \ / 
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(136) 



whose expectation value at 5 = is zero for it consists of 
cross products of n's, which make them point in different 
directions. The correlation function of two n's pointing 
along different directions is zero. 

A lengthy but straightforward calculation yields the 
susceptibility, which is a nonuniversal constant times uo/k 



fc2 



5ij 



(137) 



The superscript u differentiates this unstaggered suscep- 
tibility from the Neel susceptibility. 

The physical response to the magnetic field is propor- 
tional to the imaginary part of the retarded susceptiblity. 
Use of the fluctuation dissipation theorem leads to 
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(138) 



It exhibits a temperature broadened spinon pole and is 
thus directly sensitive to spinon creation. We see that the 
response is a universal scale invariant function of co/T 
and \k\/T. When the energy and momentum are just 
right, the system takes energy from the magnetic field 
by shooting off spinous. This behavior is illustrated in 
Figure [T2l At zero temperature the response is 



(139) 



On the other hand, it might be hard to apply a 
magnetic field with precisely the right relation between 
mometum and energy. Thus the static response, obtained 
by integrating over all frequencies, might be of more rel- 
evance. It is given by 




1 - e-l^l/^ 
1 1^1 



1^1 
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and illustrated in Figure [T3l 



Static Response to External Magnetic Field 




FIG. 13: Physical response of the system to an external, static 
magnetic field as a function of the ratio of momentum to 
temperature. 

We see that R^{k)/T is a universal function of \k\/T. 
At zero temperature we have 



R'ik) - \k\ . 



(141) 



A similar important quantity is the response to a time 
dependent but uniform field, obtained by integrating over 
all momenta 



1 
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o/T 
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. (142) 



R"{uj)/T'^ is a universal function of uj/T which is illus- 
trated in Figure [HI At zero temperature we have 



i?"(a;) ~a;2^(w) 



(143) 



We have seen that deconfined quantum critical points 
are examples of new types of quantum phase transitions 
where the standard Wilson-Ginzburg-Landau criteria 
that critical properties are controlled only by order pa- 
rameter fluctuation fails as has been recently suggested^. 
We have uncovered new properties and physics of such 
deconfined critical points as well as predicted some of 
their experimental features. 
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FIG. 14: Physical response of the system to an external, time 
dependent, uniform magnetic field as a function of the ratio 
of frequency and temperature. 
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